Semiclassical quantization of electrons in magnetic fields: the generalized Peierls 

substitution 



00 1 

o : 
o ■ 

(N ■ 
=5 

<: 

in ! 



Pierre Gosselin 1 , Hocine Boumrar 2 ' 3 and Herve Mohrbach 2 
1 Institut Fourier, UMR 5582 CNRS-UJF, UFR de Mathematiques, 
Universite Grenoble I, BP74, 38402 Saint Martin d'Heres, Cedex, France 
2 Laboratoire de Physique Moleculaire et des Collisions, ICPMB-FR CNRS 2843, 
Universite Paul Verlaine-Metz, 57078 Metz Cedex 3, France and 
3 Laboratoire de Physique et Chimie Quantique, Universite Mouloud Mammeri -BP 17, Tizi Ouzou, Alqerie 

A generalized Peierls substitution which takes into account a Berry phase term must be considered 
for the semiclassical treatment of electrons in a magnetic field. This substitution turns out to be an 
essential element for the correct determination of the semiclassical equations of motion as well as for 
the semiclassical Bohr-Sommerfeld quantization condition for energy levels. A general expression 
for the cross-sectional area is derived and used as an illustration for the calculation of energy levels 
of Bloch and Dirac electrons. 

PACS numbers: 



i 

C 

O 

o 



> 
in 
in 

(N 

in 
o 

oo 
o 



- 1—^ 

x 



Semiclassical approaches are very important in many 
area of physics for the study of the short wave length be- 
havior of quantum systems, including Bloch electrons in 
crystals or Dirac particles in external fields. An essential 
ingredient of these approaches is the Bohr-Sommerfeld 
quantization condition, whose generalization from scalar 
to vector wave fields has revealed new gauge structures 
related to Berry's phases jlj. This paper presents a de- 
tailed study of the semiclassical quantization for a sin- 
gle quantum particle in a magnetic field, exemplified by 
electrons in a crystal and by Dirac electrons. This uni- 
fied description of a particle in a magnetic field is based 
on a method of semiclassical diagonalization for an arbi- 
trary matrix valued Hamiltonian developed previously [2j 
(for a generalization to higher order in h see Q). This 
method results in an effective diagonal Hamiltonian in 
terms of gauge-covariant but noncanonical, actually non- 
commutative, coordinates, ft will be shown that a gen- 
eralized Berry's phase dependent Peierls substitution is 
necessary for the establishment of the full equations of 
motion including Berry's phase terms. This substitution 
turns out to be also an essential ingredient for the Bohr- 
Sommerfeld quantization condition of an electron in a 
magnetic field. Indeed, when reformulated in terms of 
the generalized Peierls substitution, this condition leads 
to a modification of the semiclassical quantization rules 
as well as to a generalization of the cross-sectional area 
derived independently by Roth Q and Fal'kovskii [j| in 
the context of Bloch electrons. 

Semiclassical diagonalization. Let us consider a system 
of a quantum particle in an uniform external magnetic 
field B = V x A described by an arbitrary matrix val- 
ued Hamiltonian H (II, R), where II = P + eA(R) is the 
covariant momentum and e > is the electric charge. We 
assume that the system can be separated into two con- 
tributions such that H (II, R) = H m (II) + y>(R) where 
H m (TV) is the pure magnetic part and <p(R) is the ex- 
ternal electric potential. In this paper we will be mainly 



interested by the magnetic contribution. The exact di- 
agonalization of this matrix valued operator through an 
unitary matrix U (II) is in general not known, and in this 
paper we apply a recursive diagonalization procedure de- 
veloped previously by two of the authors. This procedure 
is based on a series expansion in the Planck constant 
of the required diagonal Hamiltonian By diago- 

nal Hamiltonian it is meant a matrix representation with 
block-diagonal matrix elements associated with energy 
band subspaces. The method is based on the knowledge 
of the zero-order diagonal representation e = UqHqUq 
where Uq is the zero-order transformation matrix, Hq the 
Hamiltonian H m , in which the components n 1 are con- 
sidered formally as classical, and therefore commuting 
operators. Quantum corrections are then re- introduced 
to yield the expression for the diagonal Hamiltonian 
Hd = U (II) HU + (TV) which, if we limit ourselves to 
the semiclassical order (the semiclassical condition being 
that the radius of curvature of the orbit is large in com- 
parison with wavelength) , has diagonal operator elements 
labelled by the energy index n which reads : 



(H d ) 



(ir n ) + <P( v r, 



ehM.B 



(1) 



Here e n (it) is the zero-order matrix element of e (it can 
itself be a matrix as for a Dirac Hamiltonian, in which 
case a block-diagonalization is considered) in which clas- 
sical variables are now replaced by the quantum opera- 
tors 7r„ = II + hA^ and r„ = R + hA r , where we have 
defined the Berry connection as being the projection on 
the nth energy band A r /^ — V n [y^r/ir] °f the matrix 
A r /n = ±t [t/V 7r / r [/ + ] . ft turns out that in the mag- 
netic case the Berry connection in momentum has the 
following expression A„ = —eA r x B. Looking at Eq. (Q]) 
we have proven that instead of the Peierls substitution 
which amounts to replace the canonical momentum 
P by the covariant one II in the energy band e n , one has 
to consider a generalization of the Peierls substitution via 
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the noncanonical covariant momentum 

tt„ = n ehA x B (2) 

(where we now use the notation A = A r ). The last 
term in Eq. (Q]) is the coupling between the uniform 
magnetic field and the magnetic moment defined as 

M (tt) = §P„ ([e,A] xA) = \V n [a x A \ . Note that 

it is common (especially in solid state physics 0) to 
write the matrix elements of the components of M as 

M i = ifcV^ j2 ( - Aj) ; ti " L ( ^ fc) — (where we used A nm = 

jr (e„ — e m ) Anm ) j which thus depend on the band-to- 
band matrix element of A. 

The appearance of the Berry connection allows us 
to define naturally non-Abelian (in general) Berry cur- 
vatures Qij (7r) = d r iAj — d r jAi + \A t ,AA where for 
simplicity we omit now band indices. Position oper- 
ators then satisfy an unusual non-commutative alge- 
bra [r 1 ,^] = ift 2 y . The generalized covariant mo- 
mentum satisfy an algebra [71"% 7r J ] = —iehe % ^ k B^ + 
ie 2 h 2 e lpk e jl}l Q pq B k Bi slightly corrected with respect to 
the usual one [IT, IP] = —ieHe^ B^ by a term of order 
H 2 B 2 which can in general be neglected. The Heisenberg 
relations between the coordinate and the momentum 
[VjTT 7 ] = ih5 lJ + ieh 2 e : ' lk Q d Bk is also slightly changed 
but by a term of order h 2 B. This contribution which is a 
direct consequence of introducing the generalized covari- 
ant momentum was overlooked in previous works, with 
the exception of Bliokh's work on the specific case of the 
Dirac equation [j| . It turns out that this term is essential 
for the determination of the genuine semiclassical equa- 
tions of motion which are 

r = dE/d-K — h-k x @(tt) 
-k = -eE - ef x B (3) 

where we defined £ = e — efiM.B. As consequence of 
the non-commutative algebra, the velocity equation is 
corrected by an anomalous velocity term tt x 0, where 
the vector defined as &i — e^O^ /2 is the Berry 
curvature of an electronic state in the given nth band, 
associated to the electron motion in the nth energy band. 
These equations of motion, where first derived in solid 
state physics context in (see also [Io|) by considering 
the evolution of the wave packet of a Bloch electron in 
an electromagnetic field. In this picture, it is the mean 
over wave packets of the operator r corresponding thus 
to the wave-packet center r c and the mean of 7r giving 
the mean wave vector 7r c that are the variables in Eq ([3]) . 
The operatorial approach reveals first that the operator 
7T is in fact a generalized covariant momentum operator 
which replaces the Peierls substitution, and second, that 
the operatorial equations of motion are not restricted to 
Bloch electrons in a magnetic field but are valid for any 
physical system described by an arbitrary matrix valued 



Hamiltonian of the kind H (II, R) = H m (II) + <p(R). In 
particular they are also valid for Dirac particles moving 
in an electromagnetic field. 

Note that in solids, for crystals with simultaneous 
time-reversal and spatial inversion symmetry, the Berry 
curvature and the magnetic moment vanish identically 
throughout the Brillouin zone 0. This is the case for 
most applications in solid state physics, but there are sit- 
uations where these symmetries are not simultaneously 
present as in GaAs, where inversion symmetry is broken, 
or in ferromagnets, which break time reversal symme- 
tries. In the same way, the presence of a strong mag- 
netic field, the magnetic Bloch bands corresponding to 
the unperturbated system breaks the time inversion sym- 
metries. In all these cases the dynamical and transport 
properties must be described by the full equations of mo- 
tion given by Eq. ([3|). In the case of Dirac particles, both 
the Berry curvature and the magnetic moment are non 
zero and the full equations of motion have to be consid- 
ered 

Bohr-Sommerfeld quantization. Having shown the ne- 
cessity of the generalized Peierls substitution for the de- 
termination of the semiclassical equations of motion, we 
now investigate the relevance of this new concept at the 
level of the semiclassical quantization of the energy levels 
for an electron motion in an external uniform magnetic 
field only (<p = 0), so that Eq. J3) becomes f = D 

and tt = -eD (||xB) with D- 1 = 1-eKB®. For con- 
venience, B is chosen to point in the z-direction B =_Bk. 
Consequently the orbits satisfies the conditions £ =const 
and 7r z =const. The semiclassical quantization of energy 
levels can be done according to the Bohr-Sommerfeld 
quantization rule 

I P±dR ± = 2iTh (n + 1 /2) (4) 

where P± is the canonical momentum in the plane per- 
pendicular to the axis tt z = cte. The integration is taken 
over a period of the motion and n is a large integer. 
Now, it turns out to be convenient to choose the gauge 
A y = BX, A x = A z = 0. In this gauge, one has 
tt z = P z = cte, and the usual covariant momentum 
Uy = P y + eBX. As BX = B(x ~ KA X ) the general- 
ized covariant momentum defined as -k v = H y + ehBA x 
becomes 

TT y = Py + eBx (5) 

which is formally the same relation as the one between 
the canonical variables, but now relating the new co- 
variant generalized dynamical operators. This relation 

with the help of the equations of motion gives P y = 
TT y — eBx = thus P y is a constant of motion so that 
§ PydY = Py § dY = and Eq.Q becomes simply 
§ P x dX — 2irh(n + l/2). Now using the definition of 
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the generalized momentum P x = 7r x + 
differential of the canonical position dX 



ehAyB and the 
= dx — HdA x = 



-t^ — hdA x , the Bohr-Sommerfeld condition Eq. (0J be- 



eB 

comes 



it x dity — 2irheB 



li,/7T - 



(6) 



where the integral is now taken along a closed trajec- 



tory r in the ir space and 



2-7T ■ 



is the 



Berry phase for the orbit T. It is interesting to note 
that in terms of the usual covariant momentum (Peierls 
substitution) we have instead of Eq. ([6]) the condition 
<fH x dHy = 2iTheB (n + 1/2) . The integration in Eq. 
((6|) defines the cross-sectional area So(s,w z ) of the or- 
bit r which is the intersection of the constant energy 
surface e (n) =const and the plane tt z =const. There- 
fore the condition Eq. ([6]) implicitly determines the en- 
ergy levels e n (n z ) . Computing now the cross-sectional 
area So(£,ir z ) — Sq(s — ehM z B ,n z ) » So(e,% z ) + dS, 
with dS = § dndir± the area of the annulus between 
the energy surface e =const and the surface e + de with 

de = —ehM z B, and where dn = ^Jdn 2 + dir 2 is an ele- 
mentary length of the tt orbit. Then, as dS can be written 



dS = $ 



zhB § 



M-dK 



where the integral is 



taken over the orbit T, one has finally 

/ 1 If M z (tt) dn 



It is common to write Sq(£ , ir z ) = 2nheB (n + 7) defining 
thus the coefficient 1 ~ \ = ~<I>b ~ ^§ \dffa^ ± \ - This 
coefficient can also be written in a different form 



1 



7 2 2tt 



1 f [v X A + M] z dn 



\de/ditj 



(8) 



with v=de/dTT. Eq. © is a generalization of a previous 
result found by Roth [4[ and Fal'kovskii [5j, in the spe- 
cific context of Bloch electrons in a magnetic field. The 
connection with Berry's phase was seen later by Mik- 
itik and Sharlai [ll|- In both Q and [TTj |. the term 
[vxi + M] was written as \V n [(— + v) x A] where 

v =— + HA is the velocity operator before projection 

on a band, and II = mR, a relation valid only for a 
Hamiltonian whose kinetic energy is II 2 /2m. Therefore 
Eq. ([7]) is more general and has a broader field of appli- 
cation, as it is a general result which applies for any kind 
of single quantum particle system in a magnetic field, 
including Bloch and Dirac electrons. Importantly the 
derivation provided here is new, and it turns out to be 
the result of the generalized Peierls substitution in the 
Bohr-Sommerfeld condition. 

Bloch electron. In a crystal, the Berry gauge 
A (k) is Abelian (a scalar operator), written in terms 
of the periodic part of the Bloch wave |u„(k)) as 



A (k) =i (u„(k)| 9k |w„(k)), where k is the generalized co- 
variant pseudo momentum (k = ir/h). Application of 
Eq. ([7]) for electron trajectories in a crystal with time 
reversal and spatial inversion symmetry, where it is ex- 
pected that, both and M vanish in the Brillouin zone, 
has been studied by Mikitik and Sharlai [llj. But these 
authors also pointed out the fact that the Berry's phase 
is non zero when the electron orbit surrounds the band- 
contact line of a metal, actually cj) B — ±1/2. Conse- 
quently, 7 = in this case, instead of the previously sup- 
posed constant value 7=1/2 which is commonly used in 
describing oscillation effect in metals. As these authors 
mentioned, measurements of 7 can allows the detection 
of band contact lines. 

As a simple application of Eq. (UJ consider a crys- 
tal with time reversal and spatial inversion symmetry, 
and where the Fermi surface is an ellipsoid of revolution 
characterized by two effective masses, a transverse m± 
and a longitudinal mi one. The energy levels can eas- 

ily be deduced. Indeed £ = h I ^7 + 2m7 ) an( ^ * ne 
cross-sectional area Sq(£,K z ) is a disc of radius square 
k^ = 2mi (j^/ti 2 — so that the energy levels arc 
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h K 2 

2m t 



which actually coincide with 



the exact ones because the energy levels of an harmonic 
oscillator keep their form at large n. 

Dirac electron. Let us consider the Dirac Hamiltonian 
H = a.Tl + (3m in the presence of an uniform magnetic 
field, with a and (3 the usual (4 x 4) Dirac matrices. The 
semiclassical block-diagonalization followed by a projec- 
tion on, say, the positive energy subspace, leads to the 
(2 x 2) matrix valued energy operator £ —e — eftM.B 
where e = \/tz 2 + m 2 (c = 1) and the magnetic moment 
is given by M = ^ — j, with L = tvxA representing 
the intrinsic orbital angular momentum [8] 2]. It turns 
out that for Dirac, the magnetic moment can also be ex- 
pressed as M =e®, with the curvature vector given by 
the matrix [1] 



(tt) = - 



2e 3 



(<T.7r) 7T 



with cr the Pauli matrices. Berry's connection is defined 
as A—i (+, 7r| |+,7r) where |+,7r) is two components 
spinor of the positive energy subspace. Consider B point- 
ing in the z-direction so that tt z — P z =const, with the 
goal to compute the Landau energy levels (LEL) as an ap- 
plication of Eq. (J7J . As the cross-sectional area So(e, P z ) 



is a disc of radius square 7r^ = e 



P 2 , the applica- 



tion of Eq. dJj) consists in replacing e by £ in 7r^ so that 
we have S (£ , P z ) = n (£ 2 - m 2 — P z ), which yields the 
semiclassical quantized LEL through the relation 



1 1 f M z dn 



Now from the Berry connection A = 2e[^+m) we 
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duce the Berry's phase 4> B = — § + T + 2e(e+r. 
where r = ±1 are the eigenvalues of the Pauli matrix 
<j z . Berry's phase is the sum of a topological part — \ 

and a non-topological r ( || + 2e (f+ m ) ) one - The contri- 
bution from the magnetic moment yields ^ § g^/g^ = 

2e 2e(e+m) J a term which exactly cancels the non- 
topological contribution of 4> B , so that finally 



2fi.Be n 



1 



P, 2 



It turns out in this example that the semiclassical energy 
quantization coincides also with the exact result. It is 
usually expected that for a massless Dirac particle the 
Berry's phase takes the topological value cf) B — ±1/2, as 
a consequence of the band degeneracy at zero momentum 
[l3| . This is not the case here because the magnetic field 
lifts this degeneracy as P z is not zero in, general. But 
it turns out that the magnetic moment contribution ex- 
actly compensates for the non-topological Berry's phase 
contribution. This cancellation can be easily understood 
from the expression Eq. ([8]) for the coefficient 7. Indeed 
from the equality [v x A + M] z = [^] z + § - = § 
we deduce the expected result 7=4 + 5 = or 1. 

For a two-dimensional Dirac system it is therefore ex- 
pected that the magnetic moment for massless parti- 
cles exactly vanishes, and that the Berry's phase takes 
the topological value <f> B — ±1/2. The electron motion 
in graphene is an interesting physical situation which 
illustrates this assertion. Indeed, graphene is a two- 
dimensional carbon crystalline honeycomb structure with 
inversion symmetry so that M =0. The hexagonal Bril- 
louin zone has two distinct and degenerate Dirac points 
or valleys (labelled by r ± 1) where the conduction and 
valence bands meet and the electronic excitations behave 
like massless relativistic fermions, so that (f> B = ±1/2 and 

consequently £ n — ±^j2heB (™ + 5 + §) [HI- Therefore 
the ground state is not degenerate as there is only one 
possibility to realize it n = 0. This result explains the 
peculiar quantum Hall effect of graphene [151 ]. 

Summary. We have shown that a generalized Peierls 
substitution including a Berry phase term must be con- 
sidered for a correct semiclassical treatment of electrons 
in a magnetic field. This substitution is essential for the 
determination of the full semiclassical equations of mo- 
tion, as well as for the semiclassical Bohr-Sommerfeld 



quantization condition for energy levels. Indeed, the sub- 
stitution in the Bohr-Sommerfeld condition leads to an 
expression for the cross-sectional area which in some sort 
generalizes the formula found by Roth and Fal'kovskii in 
the context of Bloch electrons in a crystal. Application of 
this formula to Dirac electrons shows the subtle cancella- 
tion mechanism between the magnetic moment and the 
non-topological part of the Berry's phase, which yields 
the Landau energy levels. 
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